An effective Hamiltonian for phase fluctuations on a lattice: an extended XY model 
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We derive an effective Hamiltonian for phase fluctuations in an s-wave superconductor starting 
from the attractive Hubbard model on a square lattice. In contrast to the common assumption, we 
find that the effective Hamiltonian is not the usual XY model but is of an extended XY type. This 
extended feature is robust and leads to essential corrections in understanding phase fluctuations on 
a lattice. The effective coupling in the Hamiltonian varies significantly with temperature. 



PACS numbers; 74.20.De, 74.20.Fg, 74.40.+k 



The physics underlying the pseudogap phenomena in 
the cuprates|0| remains unclear and is a topic of intense 
debate. One possibility is that the pseudogap is a pre- 
cursor of the superconducting gap. This picture finds 
support in ARPES experiments which show that the 
pseudogap seems to evolve into the superconducting gap 
as the temperature (T) is lowered through the supercon- 
ducting transition temperature (Tc). In this picture, the 
Cooper pairs are preformed above Tc but without phase 
coherence^, and the superfluid stiffness is destroyed by 
the phase fluctuations. A different precursor scenario 
based on an extension of the BCS theory which induces fi- 
nite momentum pairs has also had some success [Q. Since 
the superconductivity is believed to resides in the Cu02 
planes and interlayer coupling is weak in the cuprates, 
it is conceivable that the high- Tc superconductors re- 
flect importantly their two-dimensional nature. Since no 
spontaneous breaking of continuous symmetry is allowed 
in 2D||], a reasonable expectation is that the supercon- 
ducting transition, at least, in the underdoped cuprates 
is not of the BCS kind but is rather of the Berezinskii- 
Kosterlitz-Thouless (BKT) type[|. 

It is well-known that the effective Hamiltonian for 
phase fluctuations in an s-wave superconductor is equiv- 
alent to HxY = JxY J dr{V9)'^, where 9 describes phase 
fluctuations. On a lattice, one has to either discretize 
HxY or derive the effective Hamiltonian from first prin- 
ciple. A common recipe to discretize Hxv is the as- 
sumption that the continuum limit still holds and, there- 
fore, the lattice version of the XY Hamiltonian [TLxy) 
has been uncritically accepted to describe phase fluctua- 
tions. In this paper, we derive an effective Hamiltonian 
for phase fluctuations in an s-wave superconductor from 
first principle. We find that the effective Hamiltonian is 
not of HxY type but rather it is an extended XY; how- 
ever, its critical behavior is still BKT-likc. In contrast 
to the usual phenomenological assumption, the effective 
coupling constant in the Hamiltonian now depends sig- 
nificantly on temperature. 

We begin with a 2D attractive Hubbard model on a 
square lattice as the simplest case conceivable 
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where Cia is a fermion field with spin a, fj, is a chem- 
ical potential, U{> 0) is the pairing interaction, t and 
t' describe, respectively, the nearest neighbor (n.n) and 
the next nearest neighbor (n.n.n) hopping. The symbol 
< ij > means a sum over the n.n pairs and ^ ij 3> 
indicates a sum over the n.n.n pairs. We will set the 
lattice constant a = 1 and also use units such that 
h = ks — I- The partition function (Z), base on the 
Hamiltonian 7i and written in the language of path in- 
tegral, is given hy Z = J VC'^VC exp[—S], where the 
action S — J drj^.^ C^drCia + H}. The range for the 
integral over the imaginary time t in the action is from 
to 1/T. We first consider only n.n hopping (t) with t' = 
and derive an effective local Hamiltonian for the phase of 
the order parameter. At finite temperature we obtain, in 
this case, an extended XY model Hamiltonian which in- 
cludes not only the n.n spin-spin interaction but also the 
second and third neighbor interaction. Inclusion of the 
n.n.n hopping makes this extended feature of the effec- 
tive Hamiltonian even more robust in the sense that this 
property is now manifested even at zero temperature. 

Introducing the Hubbard-Stratonovich 

transformation with an auxiliary field — Aje"^^^ , 
and making a gauge transformation '^j — exp[iT39j /2]xi, 
where Xi is a spinor for neutral fermions and are 
Pauli matrices, after integrating out the fermion fields, 
we obtain Z = J V(j)*V4> exp[— 5*6//], where an effective 
action is 
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Here, Tr means the trace over the functional space, de- 
pending on the representation, and spin space. We will 
use a symbol tr for a trace over spin space only. In 
the real (lattice) space, the Green function can be repre- 
sented as G~^{i,j) = G'^^{i,j) — , with 

Gg'^iiJ) = i-Todr + T3^J■ + T2A)Si,j +T3t^Sj^i+S (3) 

.5 



and the self energy has the form 



foSW(*,j) + ^3S(''(*,j), (4) 
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where S(o)(i, j) = it 5^- ,+5 sin(0,,,/2) and S(3)(i,j) = 
tJ2s ^j,i+s[^ — cos(0ij72)]. Here 9ij = 9i — 9j is a phase 
difference between sites i and j, and S — ±i and iy on 
a square lattice. In the expression for we consider 
only fluctuations of the phase and assume no fluctuation 
in the magnitude of the gap. Magnitude fluctuation will 
also change the self energy S. However, we ignore this 
and concentrate on the effects of phase fluctuations. For 
simplicity we consider the static case {drOi = 0). A gen- 
eralization to the time-dependent case is straightforward. 

The effective action 5'e// can be separated into a 
mean- field part (S^^^) and a phase fluctuation part 
(Se), namely, S,ff = 5^ + TrEr=i [^0^]", where 
S'(o) = JdTj2i ^1^1^ - Trln[G^i]. It can be shown|| 
that in the saddle-point approximation dS'-^yOA = 
reduces to the BCS mean-field gap equation: A — 
X^k t^'^'^ (^)' dS^-^^dfi gives an equation for 
the filling factor: n ~ 1 — J^k 'eZ ^^^^ ' ^'^^^'^ 

Ek — i/Ck + '^it^ = -2i [cos(fc:E) -I- cos(A;j^)] - fi. 
When we include the n.n.n hopping, has an additional 
term — 4t' cos(fc2:) cos(A;j,). In the calculation of the ef- 
fective phase-only action Sg, we have assumed that the 
phase fluctuations are small so that we retain only up 
to the square of the phase difference between two sites 
in expansion of So- This means that it is sufficient to 
consider only the first and the second trace for Se; in 
other words, Sg ^ S^^^ + S^^\ where S'(i) = Tr[GoI]], 
and S-t^) = iTr[GoI]GoS] . 

For simplicity, we will use four vector notation: 
K = (k,w„), = T'Ec^„Ek' where ujn is a 

fermionic Matsubara frequency, and J dXi — J dr'^^. 
It can be shown that {K\G\K') = 5{K - K')G{K') 
and {K'\'i:\K) = fot^°\K' , K) + r^t^^^K' , K) , where 
S(0)(X',if) = It ^ dX^e-'^^^'-^^^'Y.5e'^'^MX^), and 
^ t jdX,e'''^^'~^"'^^Y.sf^^'^Bs{Xi) with 
A!,{Xi) = sin{0,,,+s/2) and BsiX,) = 1 - cos(0,,,+5/2). 
Since S'^^' = trJ2K,K'{K\Go\K')(K'mK), we obtain 

5(1) =t^r^)cos(fc,) fdrY, [l-cos(0.,,/2)], (5) 

K <i]> 

where r[^^ = tr [G{K)T3e'''''^"^^] with 77 ^ 0+ being a 

convergence factor]^, which originates from G(0, r — t+). 
Note that makes no contribution to S^^\ Since 5'*^^^ 
include the phase difference only between the n.n sites 
and the first non-trivial term is 0f we obtain the usual 
XY model action; however, 0f j terms will also appear in 
5^^-'. Moreover, the phase difference between the n.n.n 
sites is also obtained from the second trace. As we will 
see later, if we include the n.n.n hopping, even in S'--^^ 
one has the phase difference between the n.n.n sites. 

Similarly, S^^^ = ^trY^^^ j^M{K)t{K,K')GiK')f:iK',K), 
where S = foS'-^-' -I- faE'^^. After some manipulations. 



one arrives at 

5(2) = -J2 A^K%'^'^°HK,K')t'^"\K\K) 

K,K' 

+ IY.^k%'^^'Hk,K')±^'\K',K), (6) 

K,K' 

where A^^^^, = tr [G{K)T,G{K')n]. The mixed 
terms do not contribute to S^^^ in our consider- 
ation for a local effective action. Since in the 
momentum space As{Xi) — ^^(Q) exp[i(3Xi], 

where Q is a four vector (q, f2,„) with a bosonic 
Matsubara frequency fi™, the first term of S*^^-* is 

A similar expression can be obtained for the second 
term. Since, as we mentioned earlier, here we are 
concerned only with the local effective action for the 
phase and wish to see if the equivalence between the 
phase-only Hamiltonian and the usual Hxy found for 
the continuum limit also holds on a lattice, we expand 
S'^^^ about Q ~ and keep the leading order. This means 
that we neglect the so-called Landau terms associated 
with damping effects, which have been shown ]lO| not 
to be important in a case of an s-wave superconductor 
in the continuum limit for T ^ 0.6Tmf, where Tmf is 
a temperature at which A (Tmf) = 0. In the s-wave 
case, the importance of the damping effects appears 
when there is a large suppression of A(T) so that 
A(r) is comparable to T, and the fraction of thermally 
excited quasiparticles is not negligible. However, it 
has been pointed out that for a d-wave superconductor 
the Landau terms have strong effects even at the low 
temperature because of the nodal structure of the order 
parameter fill . It can be shown that the first term in S*'^' 

is ¥^ J:5,5' Ek A^°^e*-(^-*') / dXas{X,)As'{X.), 
and the second term becomes 

Em' Ek trAP^e^^<'-'') J dX,Bs{X^)Bs. (X,). 
In order to compare the effective Hamiltonian with the 
lattice version Ti.xY, we need to expand S'^^^ in terms 
of the phase differences while keeping terms up to 6^ y 
Since the first non-trivial term for the second term is 
df y we will ignore it in the effective Hamiltonian. 

Using symmetry properties of G{K) such as remaining 
the same with respect to the exchange of and ky, we 
obtain the effective local Hamiltonian: Hg = {9\M\9), 
where (6*1 = [9i,i+x,di,i-x,di,i+y^(^^,i-v) and 

/ a /3 7 7 \ 

M= ^ a 7 7 
7 7 a /3 
\l 1 P a J 

with components as follow: a = ^t^^v'^^ cos{kx) + 
¥'Ek^k%' P = \t^EK^%%^os{2kx), and 7 = 
X^if ^^K^K cos(kx) cos(fcy). As one can see, the 4x4 
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matrix M is not diagonal. This means that the effec- 
tive Hamiltonian Tig is not equivalent to the usual Hxy, 
which would be pure-diagonal; namely, we have not only 
terms like Of^^- and 9'fi_^ but also terms of the form 
9i,i+x0i,i-x in He- We will next investigate effects of 
these off-diagonal terms. 

It is worthwhile making sure our new result for Tig 
reduces to the well known XY-type Hamiltonian|^ in the 
continuum limit. For this purpose, we need to recover the 
explicit value of the lattice constant a in the expression 
for Tig in order to track orders of a in the limit a —> 
and the phase difference becomes a derivative of the 
phase (V6') while ta^ where m is an effective 

mass of an electron. It is straightforward to show that, 
in the continuum limit with — ^ — fJ-i T^e ^ H^'' + 
Hf\ where H^^'^ = i^T. k'^^'k^ S dv Of , and i/f^ = 

^K%^^ I d^ (V6')^. Consequently, the effective 
Hamiltonian Hg we derive does reduce to Hxy- 

Let us now consider effects of the off-diagonal terms 
in Tig. At T = 0, Hg again becomes equivalent to the 
XY-type Hamiltonian because the off-diagonal terms of 
M, P and 7 vanish. However, at finite temperature they 
are nonzero so that Hg is no longer of the usual XY type. 
Instead 

Hg = a ^ 9f j + 2/3^ {di,i+xOi,i-x + (^i,i+ydi,i^y) 

<ij> i 

+ '^j''^{Oi,i+x + Oi^i-x) {di,i+y + di^i-y) . (7) 
i 

Since the phase fluctuation between two sites is small by 
assumption, one can make the approximation; Oi^i+x — 
sva{di_i+x)- Introducing a 2D classical spin = 
(cos(0i), sin(0i)) at site i, it can be shown, within 
the approximation we have made, that, for example, 
di,i+xOi,i-x — (Si X Si+x) ■ (Si X Si_s). Apply a vector 
identity to this expression, we then obtain (Si x Si+s) • 
(Si x Sj_;j) = SiJ^x ■ Si_j — (Si • Si+g)(Si • Si—x), where 
we used the fact that Si • Si = 1. As one can see, a 
spin at the site i -\- x couples to a spin at i — i, which 
is a n.n.n spin-spin interaction. Note that the only as- 
sumption we have made is that terms higher than Of ^ 

for a given (i,j) are negligible. This does not mean, 
however, that only the n.n phase differences are impor- 
tant. Following the procedure we briefly described, one 
can show that the second term of Hg, which is pro- 
portional to /3, becomes /3l]<ij> ^Ij + PH<ij>s. ^» ' 
within a constant, where the symbol < ij >3 indicates 
a sum over the next-next nearest pairs. Similarly, the 
third term of Hg proportional to 7 turns out to be 
27E<,j> ^Ij + 27E«y» S^ • Sj. Since Of^ can also be 
represented in term of Si • Sj, the effective Hamiltonian 
Hg can be written as 

He = -Ji ^ S, S,+,72 S,-S,+J3 ^ S,-S, , (8) 

<ii> <.i3:S> <ii>3 

where Ji = 2{a + P + 27), J2 = 27, and J3 = /3. It 
is clear that this Hamiltonian is not of the usual XY 



type but instead of an extended XY type. A geometri- 
cal explanation for the appearance of the n.n.n and the 
next-next nearest neighbor (n. n.n.n) term is illustrated 
in Fig. 1(a). Indeed, these terms come from the sec- 
ond trace proportional to t^, which has a factor E^^^S^^). 
Each self-energy E^"^ picks up (5 = ±x, ±y, and the sec- 
ond trace gives terms with resulting vectors S + S' ^ S2 
or 63, where 62 = ±x ± y and ^3 ±25:, ±2y. This ge- 
ometrical picture also works when we include the n.n.n 
hopping (t') in the electron dispersion curves. 

The physics of Hg depends on the relative magnitudes 
of coefficients Ji, J2, and J3 as well as their relative signs. 
For example, if /? is negligible, and a and 7 are both pos- 
itive so that Ji > 2J2, Hg describe a non- frustrated XY 
model, and its critical behavior can be understood 
in terms of the usual XY Hamiltonian with an effec- 
tive coupling constant Jeff = (Ji — 2J2). However, 
this does not mean that in this case Hg is equivalent to 
—JeffJ2<ij> Sj • Sj because the local behaviors of these 
two Hamiltonians are different. In general, however, as 
long as Ji is dominant, the large length scale behavior 
is of the usual XY type. To calculate Ji, J2, and J3, 
we need to know how A and fj, change with increasing 
temperature. We choose the pairing interaction U = lAt 
and the filling factor n — 0.9 and self-consistently solve 
the equations for gap and filling factor to determine A(T) 
and /i(T). What we obtain in the numerical calculation 
is that i) J3 is negligible from T = to T = Tmf, and 
ii) the inequality Ji > 2 J2 holds almost all the way up 
to T = Tmf- Near Tmf, Ji can be a little less than 
2J2; however, as we mentioned before the local effective 
Hamiltonian for phase fluctuations may not be valid in 
this temperature regime, where effects of thermally ex- 
cited quasiparticles cannot be neglected. Consequently, 
for the interesting temperature range Je// = Ji — 2J2; 
namely, Jeff — 2{a + P). Interestingly, Jeff does not 
depend on the coefhcients of the n.n.n and the n.n. n.n 
interaction terms. However, this is not the case if we 
include the n.n.n hopping. 

When we include the n.n.n hopping (i'), the Green 
function Go{i,j) and the self energy change as follows: 

G^\i,j) ^ G^H^,i) + ht'Y.sJ,,^+S2, and E(o)(z,j) ^ 
E^°) [i,j)+it' 5j^i+s2 siii(0i,j/2), and a similar expres- 
sion forE(3)(i,j) holds. In the same way that the term 
gives interactions of the form Si • Si+^j and Si • Si+53 , the 
tt' term induces the form Si ■ Si+5^ because 5 + 82 = (^4, 
and the t'^ term gives rise to Si • Si+^g and Si ■ Si+^j be- 
cause 82+52 = 5^ or (^5, where 5^ = ±a;±2y, ±2i;±y, and 
5^ = ±2x ± 2y. Geometrical descriptions for these terms 
are also presented in Fig. 1(b) and Fig. 1(c). Following 
the same manipulation as before, we obtain the effective 
Hamiltonian including the n.n.n hopping contribution: 

He = — l7i ^ ^ Si ■ Sj + ^2 ^ ^ Si ■ Sj + J-i ^ ^ Si ■ Sj 

<ij> <ij>3 

+Ja ^ S. ■ S, + J5 ^ S., • S, , (9) 

<ii>4 <ij>5 
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where ^ 



<y>4(5) 



It is obvious that the 



extended features of Tie are robust in band structure pa- 
rameters. Moreover, contrary to what we found before, 
now J2 is finite even at T = because the n.n.n hop- 
ping {t') contributes to the first trace as well as to the 
second trace||l^. In order to calculate Ji {i = 1,2,3,4) 
we choose t' = —0.2t with the same values of parameters 
as before {U — lAt and n — 0.9). Since our numeri- 
cal calculation indicates that J7i and J'2 are dominant 
and Ji > 2J2 in most of the temperature range, we ne- 
glect J3, J4, and J5. In this case as we mentioned ear- 
lier, the critical behavior is still of the usual XY type so 
that we introduce an effective coupling Jef f = Ji ~ ^Ji- 
In Fig. 2, we plot ^Jeff vs T scaled by Tmf- Since 
Jeff < near Tmf, one might think this implies a frus- 
trated XY case; however, we again point out that A(T) 
is comparable to T in this region and, therefore, we es- 
timate the validity of the local effective Hamiltonian to 
be < r < Te, where Tg is a temperature such that 
Jeff > 0. We also considered a case of a stronger inter- 
action {U = 2t) and a lower filling factor (n = 0.4) with 
the same t'. The behavior of Jeff is similar to the one in 
Fig. 2. We emphasize that Jeff approaches zero earlier 
than does the gap A(r). Note that in the continuum 
limit 1^ the coupling constant vanishes at Tmf- 

In order to find out the BKT transition tempera- 
ture {Tbkt), one has to solve self-consistently Tbkt = 
^JeffiTBKr) with equations for the gap and the filling 
factor. Since A(r) and fJ,{T) have been self-consistcntly 
obtained for given parameters, the solution of T = 
^JeffiT) (the crossing point between the dashed fine 
and the solid curve) gives Tbkt, which is indicated by 
an arrow in Fig. 2. Suppose one use the usual XY Hamil- 
tonian including only Ji term in Eq. ^, as seen (dotted 
curve) in Fig. 2, Tbkt is larger than Tmf, which is un- 
reasonable. Corrections to the usual XY model are es- 
sential in our approach. 

In summary. An effective Hamiltonian for phase ffuctu- 
ations has been derived starting from an attractive Hub- 



bard model on a lattice. Unlike the common assump- 
tion, we found that the effective Hamiltonian is not of 
the usual XY type but is extended XY. This extended 
feature is robust and has important consequences. It 
is reinforced when the next nearest neighbor hopping is 
considered. The critical behavior is still BKT-like but 
with an effective coupling which depends significantly on 
the temperature. 

This work was supported in part by the National Sci- 
ences and the Engineering Research Council of Canada 
(NSERC) and by the Canadian Institute for Advanced 
Research (CIAR). 



[1 
[2: 

[s: 
[4: 
[s: 
[6: 
[7: 

[8 
[9 

[10: 
[11 

[12 
[13 



Rep. Prog. Phys. 62, 61 (1999) 
51 (1996); A. G. 



T. Timusk and B. Statt, 
and references therein. 

H. Ding et al. , Nature (London) 382, 
Loeser et al. , Science 273, 325 (1996). 

V. J. Emery and S. A. Kivelson, Nature (London) 374, 
434 (1995). 

Q. J. Chen, I. Kosztin, B. Janko, and K. Levin, Phys. 
Rev. Lett. 81, 4708 (1998). 

N. D. Mermin and H. Wegner, Phys. Rev. Lett. 17, 1133 
(1966). 

V. L. Berezinskii, Sov. Phys. JEPT 32, 493 (1971); J. M. 

Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 (1973). 

See, for example, J. W. Negele and H. Orland, Quantum 

Many-Particle Systems (Addison-Wesley, 1988). 

V. M. Loktev, R. M. Quick, and S. Sharapov, Phys. Rep. 

349, 1 (2001) and references therein. 

See, for example, A. A. Abrikosov, L. P. Gorkov, and I. 

E. Dzyaloshinski, Methods of Quantum Field Theory in 

Statistical Physics (Dover, New York, 1975). 

I. J. R. Aitchison, G. Metikas, and D. Lee, Phys. Rev. B 
62, 6638 (2000). 

S. G. Sharapov, H. Beck, and V. M. Loktev, Phys. Rev. 

B 64, 134519 (2001). 

P. Simon, J. Phys. A 30, 2653 (1997). 

We will present detailed expressions for J7i (i = 1, 2, 3, 4) 

in a later publication. 



FIG. 1 ( Kim and Carbotte ) 



5 













52 


(a) 








■/ 













6 




(b) 








< 










5 


















(c) 








53 












► 

















FIG. 1: Geometrical illustrations of induced interactions for 
the next-nearest and the further-neightbor spin pairs. 




FIG. 2: ^Jeff (solid curve) as a function of temparature (T) 
scaled by Tmf- The BKT temperature Tbkt is indicated by 
an arrow. If one use the usual XY Hamiltonian including only 
Ji term, Tbkt is larger than Tmf (See the dotted curve). 



